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5 

TECHNICAL FIELD OF THE INVENTION 

The present invention relates to non-invasive cardiac and respiratory 
monitors. More specifically, diagnostic tools and metfiods for evaluating the 
heart and respiratory system and Improving the function of the falling heart are 
10 disclosed. 



BACKGROUND OF THE INVENTION 

Cardiopulmonary monitoring of clinically unstable patients is essential 
15 for early detection of potentially life threatening changes in the functionality of 
the heart or respiration. Cardiopulmonary monitoring is presently done on 
hospitalized patients in critical care, such as intensive care units (ICU), coronary 
care units (CCU) and perl-operatively as well as in other specialized sections of 
the hospital. Certain types of cardiopulmonary monitoring are also performed on 
20 patients outside the hospital, such as patients who suffer from asthma, high 
blood pressure, or cardiac arrhythmias. Cardiopulmonary monitoring is also 
used during exercise performance tests of patients and athletes. Clearly, 
monitoring and early detection of impending catastrophes are highly desirable in 
clinical medicine. 
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Continuous monitoring of cardiopulmonary wellbeing is presently focused on 
physiological parameters that characterize, individually, the activity of the heart 
and respiration. These parameters include monitoring of the electrocardiogram, 
blood pressure in the arterial, venous and pulmonary circulations and cardiac 

5 output for assessing the cardiovascular system. The ainn/ay pressure, 
respiratory rate, tidal volume, flow rate, pulse oximetry, exhaled CO2, and 
esophageal pressures are the respiratory parameters that are commonly 
monitored. More recently continuous overnight monitoring of breath sounds was 
introduced to detect and quantify wheezing activity in asttimatics. Additional 

10 methods were introduced to monitor patients inflicted with congestive heart 
failure (CHF) to identify early decompensation. These metiiods are based on 
continuous monitoring of the electrical impedance of the thorax, but have not 
generated reliable levels of sensitivity and specificity. Monitoring of respiratory 
crackles has also been used as an early sign of lung congestion. 

15 The activities of the heari: and the lungs are well known to be closely 
interrelated. The heart rate, blood pressure and blood flow into and out of the 
heart are influenced by the breathing cycle. The pressure inside the chest 
becomes more negative during inspiration to enable inflow of air into the iung 
alveoli. This sub-atmospheric pressure also affects the heari:, blood vessels and 

20 blood flow. In particular, the increased intrathoracic negative pressure expands 
the right atrium and ventricle, dilates and elongates tiie blood vessels, amplifies 
the ventricular filling and changes tiie position of the inter venti^icular septum. In 
addition, the changes in lung volume modify the afferent neuronal activity in the 
vagus nerve, leading to modulation of the heart rate during respiration. During 

25 quiet breathing the intiapleural pressure decreases from about -3 mm Hg to 
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about -6 mm Hg. This causes dilation of the intrathoracic segment of the vena 
cava, increased venous return to the right atrium and ventricle. The increased 
diastolic filling of the right ventricle amplifies its stroke volume into the pulmonary 
circulation by the well-known Starling mechanism. At the same time, the 
5 inspiratory displacement of the abdominal content by the contracting diaphragm 
increases the intra abdominal pressure, which further pushes blood into the 
thoracic cavity. The intrathoracic pressure changes during breathing also 
influence the output of the left ventricle, but to a lesser degree. During deep 
breathing, or in pulmonary diseases that affect the mechanical properties of the 

10 airways and the lung parenchyma these, respiratory swings of flow, resistance 
and volume are greatly exaggerated (Mountcastle V. O. Medical Physiology, 
12^ ed. Mosby Company, St Louis 1968). 

Positive pressure ventilation squeezes the pulmonary capillaries and 
increases the resistance to blood flow through the lungs. These pressures may 

15 diminish the output of the right ventricle due to increased afterload, while at the 
same time reducing the output of the left ventricle due to fall in its diastolic filling 
and preload. These phenomena are known to result in wide fluctuations of the 
cardiac output and blood pressure. Similar fluctuations are also seen during 
resisted breathing, cough and isometric muscle straining. 

20 Monitoring of heart sounds is well known to medicine and physiology 

for many years, even before the invention of the stethoscope by Laennec in 
1819. The first and second heart sounds are associated with closing of the atrio- 
ventricular and ventricular outlet valves, respectively. They are loud and distinct 
sounds with somewhat different amplitude and template in different areas of 

25 precordial auscultation. The first and second heart sounds are mcxlified, and 
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sometimes are even completely missing during diseases of the heart or the 
lungs. Additionally, the 3^^ and 4^^ heart sounds are well known to be 
associated with defects of the left ventricular filling during diastole. Rumbling 
sounds heard in between the first and second sounds (systole) or the second 

5 and first sounds (diastole) are appropriately called systolic and diastolic 
mumiurs, respectively and are associated with abnomnal blood flow through the 
nan-owed or malfunctioning heari: valves. Infomiation on heari: sounds, 
phonocardiography and the art of interpretation of heart sounds is readily 
available in many text books and articles, such as in "Rapid interpretation of 

10 heart sounds and munnurs" by Emanuel Stein, Abner J. Delman (Editors), 
Williams & Wilkins, 1997, 4th edition, the contents of which are incorporated 
herein by reference. 

The respiratory changes of cardiac activity are well known. In 
particular, the phenomenon of "Pulsus Paradoxus" is a recognized sign of 

15 severe asthma and ainA^y narrowing. It is defined as a decline of greater than 
12 mm Hg (in some texts 20 mm Hg) in the systolic blood pressure during 
inspiration. Detection of Pulsus Paradoxus in a dyspneic patient is an ominous 
sign that calls for aggressive and immediate intervention. On the other hand, 
cardiac anhythmia, the accelei^tion and deceleration of the heart rate during the 

20 respiratory cycle is often a benign and normal phenomenon, particularly in 
young children. The only available infonnation on the effect of respiration on the 
heart sounds per-se is on the width of splitting of the second heart sound and on 
certain cardiac murmurs. OthenAnse, no infonnation is available on the changes 
in the heart sounds induced by respiration and on the extent of these changes 

25 during various cardiopulmonary conditions. 
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The failing heart condition is a cause for severe morbidity and mortality 
rates. This condition can be helped by a broad range of assist methods. These 
include phamiaceutical agents (e.g., digitalis and other positive inotropic 
agents), mechanical support (e.g., intra aortic counter pulsation, and 

5 intravascular coronary artery stent), and various types of electrical stimulation 
(e.g., ventricular resynchronization therapy, Guidant*s Contak CD, Medtronic*s 
multi electrode epicardial pacing. Impulse Dynamics' timed refractory period 
current). All cardiac assist methods must be monitored as regards their 
performance. This feedback can detemiine if, and to v^at extent the assistance 

10 is effective. Such monitoring is either perfomied in intervals (e.g., by periodic 
detemiination of the left ventricular ejection function), or continuous (e.g., by 
monitoring the left ventricular pressure - Remon Medical Technologies LTD of 7 
Haiamish Street, Caesaria Industrial Parte, 38900, Israel). To be effective, the 
mechanical and electrical assist methods should be synchronized with the 

15 cardiac cycle so that they can be activated at specific phases of the cardiac 
cycle. For example, in US patent 6,285,906 the contents of which are 
incorporated herein by reference, impulse dynamics provides current stimuli at 
certain precise phases of the heart cycle. To do so, most methods use the EKG 
(electrocardiographic) signal to find the appropriate timing within the cardiac 

20 cycle for activation of the assistance. However, the EKG may not provide 
suffidentiy detailed or accurate information about the timing cardiac cycle, or 
may be contaminated by electromagnetic noise. 



wo 2004/035137 



PCT/IL2003/000858 



6 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1. is a schematic description of the signal collection and processing device 
providing synchronized features of the cardiac activity; 

5 Figure 2. is a schematic description of the entire process in accordance with a 
preferred embodiment of the invention, perfomiing propbabilistic analysis. 

Figure 3. is a schematic description of the decision loop involving cardiac assist 
control in accordance with the invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

The present invention provides a novel method for monitoring the intenelated 
functionality of the heart and the respiratory system by tracking the changes in 

5 cardiac sounds with respect to the activity of the heart during the respiratory 
cycle. A novelty of the new approach is that the timing of each heart contraction 
within the breathing cycle is used as a cardinal parameter in the analysis of the 
heart sounds relative to the electrocardiographic signal. By doing so, the heart 
sounds generated during similar segments of the respiratory cycle, for example, 

10 the first part of inspiration, are analysed together to generate respiratory 
synchronized stable templates and features of the heart sounds and their 
interaction v\rith the EKG. Thus, the extent of changes in EKG - modulated heart 
sounds parameters during each phase of the respiratory cycle may be tracked 
over time to detect variability of the breathing related swings that indicate 

15 worsening of the status of the cardiopulmonary system. In addition, the invention 
facilitates tracking of the stable synchronized cardiac sounds features over time 
that are indicative by themselves of the well - being of the mechanics of the 
heart and lungs. It is further disclosed that the detemiination of the stable 
synchronized cardiac sounds features is classified by the timing and template of 

20 the EKG. Together, these new methods provide valuable infomiation on the 
status and mechanics of the complete cardiopulmonary system. 

The invention overcomes two major hindrances to extraction of stable 
features from heart sounds: the variability of the cardiac sounds during the 
respiratory cycle and changes in the sounds caused by variability of the 

25 electrical activity of the heart The new method extracts features such as 
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amplitude, duration, frequency content, template and chirp components of the 
heart sounds and clusters these features with respect to their timing in the 
respiratory cycle and /or their underlying EKG morphology. This method 
presences the attributes of the cardiac acoustic signal that are otherwise 
5 eliminated or greatly diminished by grouping the sounds or their features without 
respect to their respiratory timing and EKG, as is cun-entiy learned by the 
existing art. The extraction of synchronized stable features greatly diminishes 
their stochastic variability, so that any changes in the detemilnistic components 
of the heart sounds, such as those caused by changes in the mechanical activity 

10 of the heart, are more readily and accurately detected and characterized. The 
acoustic properties of the heart sounds reflect the mechanics of the cardiac 
contraction, the integrity of the cardiac structure and dynamics of flows, volumes 
and pressures in the cardiovascular system. Therefore, detection of changes in 
the synchronised stable sounds relative to a baseline may be used to alert the 

15 patient, his or her caregivers, or an automated algorithm or device that can be 
prompted to provide a remedy to the patient's condition. This new method, while 
preserving complete non-invasiveness, provides far more comprehensive 
monitoring of the patient than the cunrent existing methods such as the EKG. 

The grouping of synchronised stable features of the heart sounds with 

20 or without respect to the EKG template defines a baseline extent of Within 
Breath Variability (WBV). Any changes of the extent of variability of the sounds 
or their features within the respiratory cycle relative to a previous period, such as 
the baseline, can be readily and accurately detected. These alterations of the 
WBV may be caused by increased pressure swings within the chest due to 

25 altered breathing pattern or changes in the resistance, compliance or other 
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mechanical properties of the lung. Such changes are often seen in asthma 
patients, patients Inflicted by congestive heart failure (Ci-IF). patients who are 
artificially ventilated, anaesthetised, or otherwise critically ill. 

5 Analysis of the heart sounds 

In Fig. 1 to which reference is now made the main stages in collecting 
and processing of the signals, in a prefenBd embodiment of the invention are 
shown. Cardiac acoustic signals 20, respiratory acoustic signal 22 and EKG 
electric signal 24 are collected by the appropriate detectors. The cardiac 

10 acoustic signal and the respiratory acoustic signal are separated and 
conditioned, by module 28, which includes in a prefenred embodiment of the 
invention filtering, denoising, amplification, multi-sensor analysis for removal of 
background noise and recording artifacts. Feature extraction module 30 of the 
invention perfonns segmentation. Heart sounds are segmented into S1, S2, S3, 

15 S4, the breathing activity is segmented into the Heart sounds are segmented 
into S1, S2, S3, S4, the breathing activity is segmented into five respiratory 
stages; first half of inspiration, second half of inspiration, first half of active 
expiration, second half of active expiration and the post expiratory pause if 
present, and the EKG signal is segmented into Its components; P, QRS, T 

20 waves. The segmented input data undergo timing analysis by module 32. 
Temporal segmentation feature parameters are extracted from the segmented 
data by module 34. The parameters are estimated individually from heart 
sounds appearing at the different respiratory stages. The EKG signal is also 
analyzed with respect to the respiratory stages. 
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Heart sounds are aoquired via a single microphone or a plurality of 
microphones located on the chest, back and/or neck of the patient. After signal 
conditioning that includes pre-amplification, the signal is band-pass filtered so 
that the respiratory component may be separated from the heart sounds. The 

5 isolated sound structures are then detected in the heart sounds signal and 
classified as SI, the first heart sound, and S2 the second heart sound. The 
degree of isolation (which is related to the beating rate) is estimated earlier from 
the data using a density plot of the inter-beat duration. A temporal window, 
related to the time interval variability, is used to distinguish between S2 and S1. 

10 Further validation of the S1/S2 segmentation may be provided by the EKG 
signal where the timing of S1 is known to correlate with the QRS sequence and 
the timing of S2 corresponds to the T wave. After the different sounds are 
segmented, an attribute conresponding to the respective respiratory cycle is 
attached to each of them, so that they can be further collected into groups with 

15 their conesponding timing within the respiratory cycle. Averaging the different 
sounds and/or their attributes in the corresponding respiratory cycle fomis 
adaptive templates. Averaging is perfomied within each timing interval on the 
raw signal, the frequency content representation, and on a time-frequency 
representation obtained by short time Fourier transform, Wigner distribution 

20 analysis or a more adapted time-frequency representation using continuous 
wavelet transfomi or best basis and discrete wavelet representation. The latter 
methods include internal denoising (compression) which is achieved by 
removing the low energy coefficients from the transfomied signal representation 
before perfomiing the inverse transform. 
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The time-frequency representation can be used for parameter 
estimation of the chirp properties (mainly S2). These properties include the 
chirp slope, namely the frequency modulation of the signal, start and stop 
frequency as well as slope, and the amplitude changes in the signal. Properties 

5 related to the time-frequency representation are extracted using conventional 
curve fitting techniques. Features related to the amplitude changes, e.g. 
increase in negative amplitude of S2, frequency content as is recorded by 
Fourier transform, and instantaneous frequency changes which is recorded in 
the time-frequency representation. The parameters corresponding to each of the 

10 above categories, and their variability may be empirically estimated from the 
data and compared to a data base or baseline values of these parameters. The 
variability together with the estimation of the mean, may further be used for the 
estimation of statistically significant changes in each of the parameters, both 
during a respiratory cycle and when compared with a nomial baseline of the 

15 specific patient or of previously recorded and labeled data. 

The separation of the detenninistic and stochastic signal components 
is done during entropy-based denoising, e.g. with wavelet or best basis 
representation. The wavelet coefficients with smallest variance, which 
conrespond to smallest contribution to the energy of the collection of heart 

20 sounds, are set to zero. This parameter shrinkage removes the part of the signal 
that is independent of the heart sounds and can be considered noise. 

A non exhaustive list of temporal segmentation parameters and 
features, with respec* to the respiratory stages extracted from the heart sounds, 
respiratory cycle and EKG useful in the cardiac assessment of the invention is 

25 provided next: 
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1. Amplitude and energy of S1 and S2 in each segment, (peak 
and RMS) where S1 and S2 are the first and second heart 
sounds, respectively. 

2. Timing of S1 and S2 relative to the EKG and relative to each 
other whereas the EKG represents the electrical activity of the 
heart, measured with skin-surface electrodes. 

3. Delta time between occurrence of S2 and S1 and occurrence of 
S1 and the S2 of the previous heart cycle. 

4. Delta time between the occurrence of S1 and occurrence of the 
QRS wave. 

5. Delta time between the occurrence of S2 and occurrence of the 
P wave. 

6. Delta time beween the occun^ence of the P wave and the R 
wave. 

7. Frequency content and instantaneous frequency content of S1 
and S2. 

8. Cord length of S1 and S2. 

9. Polarity of first and largest component of SI and S2. 

10. Amplitude incline and damping (decaying) rate constant of S1 
and S2. 

1 1 . Chirp parameters of SI and S2. 
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12. Duration of whole cardiac cycle (reciprocal of the heart rate) as 
well as of each heart sound segment (i.e. of S1, S2, S3, S4, or 
mumiurs). 

13. Signature (or shape) specific for each sensor by using signal 

s processing methods and/or polynomial fit. Differentiate between 

deterministic elements (e.g. S1^) and stochastic elements such 
as munnurs and/ or artefacts. Determine separately for each 
frequency range. Obtain template pattern of the heart sounds 
signature for each segment of the respiratory cycle. 

10 14. Cross correlation between various precordial sites. An 

expansion of this is creating a dynamic acoustic map of the 
anterior chest using many sensors, each displayed on an array 
with brightness indicating amplitude and color(s) indicating the 
frequency content. The dynamics of the changing map can be 

15 tracked for changes during the respiratory cycle. 

15. Calculate the deconvolution of the heart sounds/ ballistographic 
signal in order to estimate the mechanical activity at the source 
that generates the signs. 

16. An AM/FM Modulation decomposition of the heart sound to 
20 analyse independently amplitude changes and frequency 

changes in the heart sounds. 

17. Amplitude multiplied by the width of each heari: sound which is 
an approximation to the temporal length of the heart sounds 
times the peak amplitude of each heart sound. Energy of some 
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components of the heart sound as is measured in a 
time/frequency representation which can be adapted to the 
signals. 

Using the above parameters and features, detemriinistic and stochastic 
5 components of the features and parameters are extracted. This defines the 
acceptable boundaries for each feature and parameter. Such boundaries are 
determined by the extent to which the parameters vary whilst the patient is in 
stable condition. In a prefenred embodiment, these boundaries are defined 
by the mean and n x SD of each parameter and feature for each the 
10 respiratory stages, where n is 2 for example. In particular, the (potentially non- 
linear) boundaries of the hyper-parameter space of the full set of the 
parameters are detemiined. Thus, for example, it is possible that a patient can 
be at a point where each of the parameters and features is within the defined 
boundaries, but the full set of parameters is outside of the boundaries. 

15 A probabilistic model, which is estimated in a prefenred embodiment, 

via calculation of maximum likelihood, is applied to the processed data to further 
refine the model estimation. This binds together the temporal progress from one 
respiratory stage to another (via a hidden Markov model). At this stage, 
additional infomnation from other sources such as blood pressure, oxygen 

20 saturation, patient weight, chest electrical impedance and other means, is fused. 
The full temporal model and set of segmented features are then analyzed to 
determine scoring of deviation firom nonnality, or from the acceptable patients' 
boundaries. 
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The probabilistic analysis assumes a Gaussian distribution to the 
continuous features of the model and a binomial (sometimes multinomial) 
distiibution to the categorical parameters. Following the estimation of the mean 
and standard deviation, significance or abnomiality of an observed features may 
5 be determined and using a probability model which takes into account the 
probability of a collection of features (taken from an already collected and soon 
to be collected data), one can estimate the probability of various abnomnal 
events. 

In Fig. 2 to which reference is now made, tiie whole sequence of steps 
10 from tiie signal transduction to model display is described schematically. In step 
50 signals are transduced to produce a coneponding electric cun^ent In step 52 
the signals are conditioned. In step 54 tfie signals are segmented to express the 
segments of tiie physiological riiytmicity. In step 56 Parameters are extracted, 
and In step 58 the detemninistic components are separated from the stochastic 
15 components. In step 60 the probabilistic analysis Is performed, accepting 
physiological data at step 62, and in step 64 scoring and display are provided. 
Altiiough not shown, the probabilistic analysis performed reflects on tiie signal 
segmentation to influence the segmentation process. 

20 Sensors for collecting signals of the heart and respiratofy system 

In general, apart from the ubiquitous EKG which senses electrical 
signals of the body, all other sensors taking part in the signal collection in 
accordance with the present invention are acoustic or mechano - acoustic 
sensors. Such sensors transduce mechanical pressure, or acoustic energy to 

25 electrical energy furt;her to be analyzed electronically. Acoustic sensors for 
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sensing activity of the heart to determine cardiac mechanical performance, 
which include typically rhythmicity, fine cardiac temporal cycle variability and 
muscle contractility, can be patched on the patients chest or can be implanted 
subcutaneously. Conditioning respiratory activity sounds and respiratory activity 
5 mainly for the detemiination of respiratory rhythm is preferably performed using 
impedance plethysmography. 

Electronic hardware for processing the data and providing 
diagnostic aids and control signals 

10 Transduced electric signals from the sensors, as well as EKG signals, 

are collected by data collecting hardware, sampled and digitized on a dedicated 
processor or in a computer. All processing and statistical calculations and model 
assessment are done in the computerized environment. A GUI or a printer is 
typically used to produce a report to the medical team in charge. In that case 

15 that control signals are required, this can be supplied by a controller linlced to 
the processor or computer wherein the processing is performed. 
IVIedical uses of the invention 

A simple example of the method of the invention includes segmenting 
the breathing cycle into exhaling and inhaling segments, and recording 

20 contracting heart sound at these inhaling and the exhaling segments. Sound 
amplitude samples are collected and then averaged, average A for the set of 
samples representing the exhaling segment period and average B for the set of 
samples representing the inhaling segment period. Statistical difference is 
calculated between the averages A and B, providing evidence regarding 

25 quantitative difference in magnitude between the two sets. The physiological 
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significance of the difference is such that the larger difference Is indicative of a 
pronounced tendency of the heart to be limited by the expanding lungs. 

These parameters and changes may be used to detect early changes 
associated with pulmonary and/or hemodynamic alterations. These include CHF 
5 (changes in the heart's mechanical activity, the lungs specific gravity (density) 
and increased fluctuations of intrathoracic pressure); mechanical ventilation 
(changes in PEEP, airway resistance), changes in hydration and blood volume 
(shock) (cardiac mechanical activity, increased breathing activity), hypertension 
(increased cardiac mechanical activity), myocardial infarction (changes in 
10 cardiac contraction sequence). 

The detection of each of the above-listed conditions or worsening 
herewith is done by comparing the respiratory synchronized cardiac sounds or 
features thereof to a patient-specific baseline, or to a template consisting of data 
from an aggregate of nomnal and abnormal patients. 

15 

Heart assist methodology 

Assessment of the mechanical condition of the heart performed in 
accordance with the invention can be used to help treating the failing heart. An 
acoustic cardiac assist method comprises, in accordance with the present 
20 invention, a detector of heart sounds consisting of a single or a plurality of 
sensors that can pick up and transduce the heart sounds of a patient into 
electrical signal and means for amplifying the signal, a detector of breathing 
activity and means for synchronizing the two signals. The heart sounds timing 
signal can be used to trigger the activity of a cardiac assist device, either directly 
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or with a delay. Such delay may be pre-determined or modified automatically, 
based on information on the heart sounds amplitude, the interval between 
them, or their amplitude and frequency content. 

in Fig. 3, to which reference is now made, an automated control 
5 mechanism a cardiac assist device using a synchronization algorithm of the 
present invention, is shown. Acoustic cardiac signals are transduced 80 and 
respiratory signals 82 are transduced to be collected by the monitor 84. 
Subsequently, they are processed at step 86 by the method of the invention to 
extract accurate parameters of the heart. If the heart is regarded failing at 
10 diagnostic step 88, the cardiac assist device is activated. 

The heart sounds timing signal Is used to actuate a demand-type or 
synchronized cardiac pacemaker. This pacemaker is chronically implanted, or 
temporarily used to defibrillate a fibrillating heart, such as used in automatic or 
manually triggered electric defibrillators. 

15 

Monitoring variability changes 

The grouping of synchronized stable features of the heart sounds with 
or without respect to the EKG template defines a baseline extent of Within 
Breath Variability (WBV). Any changes of the extent of variability of the sounds 

20 or their features within the respiratory cycle relative to a previous period, such as 
the baseline, can be readily and accurately detected. These alterations of the 
WBV may be caused by increased pressure swings within the chest due to 
altered breathing pattem or changes in the resistance, compliance or other 
mechanical properties of the lung. Such changes are often seen in asthma 

25 patients, patients inflicted by congestive heart failure (CHF), patients who are 
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artificially ventilated, anesthetized, or othenwise critically ill. When such changes 
are found they can be used to alert the patient, his or her caregivers, or an 
automated algorithm or device that can be prompted to provide a remedy to the 
patient's condition. 

5 

Short temi duration analysis 

Known in the art are changes in the duration of each heartbeat in 
association with the breathing cycle. These changes are usually detected by 
measuring the distance between consecutive "R" waves in the EKG that are 

10 called "the R-R interval". The changes in the duration of the R-R interval over 
time are called "Heart Rate Variability" (HRV). HRV is considered as indicative 
of the activity of the sympathetic and parasympathetic control of the heart. It is' 
substantially reduced In gravely ill patients. Recently, changes (reduction) in 
HRV have been associated with particulate matter (PM) air pollution and are 

IS believed to be an important mechanism in the increased mortality due to PM 
pollution. Other conditions, such as the sleep apnea syndrome, chronic cardiac 
ailments and the use of certain medications are also known to be associated 
with changes in HRV. 

In accordance with the present invention, analysis of the short-temi 

20 duration changes of systolic time (Tsys) and diastolic time (Tdias) is perfomned 
and studied with the respiratory cycle as well as the other amplitude and signal 
morphology features that have been discussed above. Of particular interest are 
conrelation between short temi changes which are indicative of a well-heart 
being such as increase in Tsys which leads to a an increase in the amplitude 

25 and sharpness of S2, indicating an increase in duration that is due to more blood 
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flow into the heart, or a negative cx)rrelation between Tsys and Tdias, indicating 
more elaborated breathing (possible air-way obstruction) but healthy heart 
dynamics, vs. changes in Tsys that are not conelated with Tdias and vice versa. 
The method of the invention can be used to differentiate between 
5 inhalation and exhalation phases of the respiratory cycle. It can also be used to 
detenmine the template and variations of S3, S4 and cardiac mumnurs. 



